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ABSTRACT 

Using high time cadence images from the STEREO EUVI, CORl and C0R2 
instruments, we derived detailed kinematics of the main acceleration stage for a 
sample of 95 CMEs in comparison with associated flares and filament eruptions. 
We found that CMEs associated with flares reveal on average significantly higher 
peak accelerations and lower acceleration phase durations, initiation heights and 
heights, at which they reach their peak velocities and peak accelerations. This 
means that CMEs that are associated with flares are characterized by higher and 
more impulsive accelerations and originate from lower in the corona where the 
magnetic field is stronger. For CMEs that are associated with filament eruptions 
we found only for the CME peak acceleration significantly lower values than for 
events which were not associated with filament eruptions. The flare rise time 
was found to be positively correlated with the CME acceleration duration, and 
negatively correlated with the CME peak acceleration. For the majority of the 
events the CME acceleration starts before the flare onset (for 75% of the events) 
and the CME accleration ends after the SXR peak time (for 77% of the events). 
In ~60% of the events, the time difference between the peak time of the flare 
SXR flux derivative and the peak time of the CME acceleration is smaller than 
±5 min, which hints at a feedback relationship between the CME acceleration 
and the energy release in the associated flare due to magnetic reconnection. 

Subject headings: Methods: statistical. Sun: coronal mass ejections (CMEs), Sun: 
flares 
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Introduction 



Solar flares and coronal mass ejections (CMEs) are the two most energetic phenomena 
in our solar system. Solar flares are abrupt releases of energy up to 10^^ J within tens of 
minutes and c an be observed in th e whole electromagnetic spectrum from radio emission 

ch 



to 7-rays (e.g. 



Fletcher et al. 



2OIII ). CMEs are sp oradic elections of coronal material wit 



Yashiro et al. 



2004 : 



Gopalswamy et al. 



veloc ities in the range of ~100-3000 km s~^ (e.g. 
20091 ). It is generally accepted that both CMEs and flares are different manifestations 
caused by magnetic reconnection in the solar corona, but the details how both phenomena 
are related are still under investigation. 

Various statistical studies using white light coronagraphic observa tions showec 
positive correlations between the flare intensity a nd CME velocity (e.g. 



Moon et al. 



Vrsna. 



2004 



c et al. 



2005 



Yashiro &: Gopalswamv 



Yashiro fc Gopalswamy 



20091 ) or CME kinetic energy ( iBurkepile et al. 



2002 



20091 ). The 



temporal differe nces between the CME and t 



f lare o nset were found to be quite small. 



ne 



Michalekl (120091 ) and 



Yashiro fc Gopalswamy 



( 120091 ) found a Gaus sian distri b ution for the difference between the flare and CME 



onsets. According to iMichalekl ( l2009l ) the mean difference is 7 min. In these studies, 
linear back extrapolation of the CME height time curve was used to estimate the CME 
start time. However, the CME onset times cannot be accurately determined by using 
only coronagraphic observations which miss the early phase of the CME acceleration and 
propagation due to the occulter disk. 



The CME kinematics typically shows three phases of evolution ( IZhang et al.ll200ll ). In 
the initiation phase, the CME rises with velocities of several tens of km s^^, followed by 
an impulsive acceleration. Thereafter, the CME propagates with almost constant or slowly 
decr easing/increasing velocity depending on its interact ion with the ambient solar wind 



(e.g. 



Gopalswamy et al. 



20001). Recent case studies (e.g. iGallagher et al 



2003 



Vrsnak et al. 
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2004J : iTemmer et al.ll2008[ ) combined EUV images with coronagraphic observations to derive 
detailed CME acceleration profiles from the CME initiation close to the solar surface until 
its propagation beyond 15 -R0. In these events the CME acceleration profiles showed a good 
synchronization with t he energy release in the associated flare, as evidenced in the HXR 



flux or SXR derivative. 



Maricic et al 



(I2OOJ) studied a set of 22 CME events associated with 
flares using SOHO EIT, MLSO Mark IV, LASCO C2 and C3. They report correlations 
between various CME and flare parameters as well as a close synchronization between the 
CME acceleration and the flare SXR flux derivative for > 50% of the events under study. 



In 



Bein et al. 



(12011 



to which we refer to as paper I in the following), we presented 
statistics and correlations between various decisive CME parameters for a sample of 95 
events: peak velocity, peak acceleration, acceleration duration, height at peak velocity, 
height at peak acceleration and initiation hei ght. To this aim we combined EUV images 



from the Extreme Ultraviolet Imager (EUVI; 



Wuelser et al. 



20041 ) with CORl and C0R2 



coronagraphic observations onboard the STEREO mission. The high time cadence and 
the overlapping fleld of views (FOV) of the different STEREO instruments enabled us to 
derive detailed and continuous CME height-, velocity- and acceleratio n-time profl l es fro m 



Bein et al. 



(120 iih 



the low corona up to about 15 R©. Out of the 95 CMEs presented in 
70 events could be associated with a GOES (Geostationary Operational Environmental 
Satellite) flare and 24 events with a fllament eruption. 9 events are associated with both, 
a flare and a fllament. In the present paper, we perform a statistical study on the relation 
between characteristic CME and flare parameters as well as on the temporal relationship 
between the two phenomena. We also study the characteristic CME parameters separately 
for events with/without associated flares and with/without associated fllament eruptions. 
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Data 



The CME kinem atics was derived from data of the STEREO/SECCHI instrument suite 



( iHoward et al. 



20081 ). The STEREO mission consists of two nearly identical spacecraft, 
STEREO- A, which moves ahead the Earth, and STEREO-B, which moves behind the 
Earth. The SECCHI package consists of an Extreme Ultraviolet Imager (EUVI), two 
coronagraphs (CORl and C0R2) and two Heliospheric Imagers (HIl and III2). For our 
study of the CME kinematics we combined data from EUVI, CORl and C0R2. 

EUVI provides observations from the solar chromosphere and lo w corona in four 



different wavelengths with a FOV up to 1.7 Rq (iWuelser et al.l 120041 ). For our study we 
used primarily 171 A observations with a time cadence as high as 75 sec. In some cases 
it was only possible to track the CME in 195 A observations with a default cadence of 
10 min but in some cases as good as 2.5 min. The high time cadence was necessary to 
derive detailed acceleration profiles of the CME initation and impulsive acceleration phase. 
Observations in 171 A and 195 A were compared to 304 A filtergrams to check if a filament 
eruption and/or flare was observed in association with the CME under study. 

The inner coronagraph CORl has a FOV of 1.4 to 4 Rq. For most events it has a 
nominal time cadence of 5 min but in some cases it was also lower (20 min). The outer 
coronagraph C0R2 has a FOV of 2.5 to 15 Rq and an observi ng cadence of 30 mi n for 



20081 ). The 



the polarized brightness images, which were used in our study (iHoward et al.l 
overlapping FOVs between EUVI and CORl and between CORl and C0R2 made it 
possible to connect the same CME structure in all three instruments, and thus to study the 
CME evolution continuously from its initation close to the solar surface up to about 15 Rq. 
To compare the CME kinematics with the flare SXR evolution we used GOES 1-8 A flux 
measurements. 
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Results 



We studied 95 CMEs, which occurred between January 2007 and May 2010, i.e. at the 
transition from solar cycle no. 23 to no. 24. For each event the full CME height, velocity 
and acceleration profiles could be derived. For details on the methods applied we refer to 
paper I. For each event, which was associated with a GOES flare, we compared the CME 
kinematics with the GOES SXR flux evolution and its derivative. Because of the temporal 
correlation between the HXR emission caused by flare-accelerated electron bea ms and th e 



derivative of the SXR emis s ion observed in ma i iy flar e events (Neupert effect 



1968 



Dennis &: Zarro 



1993 



Veronig et al. 



2002 



Neupert 



20051 ) we use the derivative of the GOES 



1-8 A flux as a proxy for the evolution of the flare energy release. Figures [T] and [2] show 
four representative events of our sample. The top panels show the CME height, the second 
panels the CME velocity and the third panels the CME acceleration curve against time. In 
the bottom panels the GOES 1-8 A flux and its derivative is plotted. 

For each CME we derived the following characteristic parameters: (1) peak velocity 
Vmax, (2) peak acceleration a^ax, (3) acceleration phase duration tacc, (4) flrst measured 
height Hq, (5) height at peak velocity h^max, (6) height at peak acceleration hamax- We 
applied a spline flt to the CME height-time proflle and from its flrst and second derivative 
we derived Vmax and amax, respectively, amax is given by the maximum value of the CME 
acceleration proflle; Vmax is deflned as the value of the velocity proflle where the acceleration 
proflle has decreased to 10% of its peak value. This deflnition is use d in order to ensur e 



that the peak velocity during the CME impulsive acceleration phase ( iZhang et al. 



2004h is 



calculated excluding the effect of the subsequent residual acceleration phase (for details see 
paper I). The acceleration phase duration tacc was deflned as the time difference between 
taccend, the time when the acceleration proflle has decreased to 10% of its maximum value, 
and taccstart, the time in the increasing phase of the acceleration when the proflle has 
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reached 10% of its maximum value, ho, the first measured height in our height-time plots, 
is used as a rough estimate of the CME intiation height, because we measured the CMEs 
as soon as it was possible to observe them in the low corona. The heights h^max and 
hamax are defined as the heights where the CME reached its peak velocity (vmax) and peak 
acceleration {umax)- A detailed study of these CME parameters is presented in paper I. 



3.1. CME characteristics in dependence of fiare/filament association 



We studied each CME parameter {vmax, cimax, tacc, Kmax and hamax) separately 
for CMEs with/without associated fiares and with/without associated filament eruptions. 
Figures [3HH] show the distributions of the different CME parameters for the whole sample 
of 95 events (grey histograms) together with the distribution for events associated with 
flares (colored histogram in the top panel) and for events associated with fllament eruptions 
(colored histogram in the bottom panel). For each parameter we derived the arithmetic 
mean together with its standard deviation and the median and its absolute deviation, for 
the whole sample as well as separately for the samples with/without associated flares or 
fllaments. The resulting values are summarized in Table [H Because the distributions do not 
show a Gaussian behaviour but have a tail towards high values, a lognormal distribution 



bette r describes their behaviour ( iLimbert et al. 



2001 



Yurchyshyn et al. 



2005 



Bein et al. 



201 ll ). Thus we also calculated the mean value /i and the stan dard deviatio n a of the 



natural logarithm of the different CME parameter distributions ( j Cowan 



19981 ): 



a 



In ( + 1 



fi = ln{E) - - (2) 

with E the mean value and Var the variance of the quantity. In the following we use 
/i* = and a* = as the median and multiplicative standard deviation of the lognormal 
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probability function of the distribution (iLimbert et al.ll200l[ ). The confidence interval of 
68.3% is given by e'^^"; /u and a are also summarized for each CME parameter in Table [H 

To check if the CME distributions with/without associated flares and with/without 
associated filament eruptions, respectively, come from the s ame continuous distribu tion 



Young 



1977 



Sachs 



19971) • This 



(null hypothesis) we used the Kolmogorov-Smirnov test (e.g. 
test is nonparametric, i.e. it is not necessarily required to know the type of distribution 
(e.g. gaussian, lognormal, exponential). A test statistic D is calculated by 

Fi F2 



D = max 



(3) 



rii n2 

with Fi and F2 the cumulative distribution functions of both subsamples and rii and n2 the 
number of events in each sample. D is then compared with 



Dr. = Kn 



111 + 77,2 



rii ■ 77,2 

where Kq, can be found in tables for different significance levels a. a represents the 
probability that we wrongly reject the true hypothesis or accept the false hypothesis 



(4) 



( lEssenwanger 



19761 ). If D>Da the null hypothesis can be rejected, i.e. both subsamples do 



not come from the same population. 



3.1.1. CME - Flare Relation 

For 70 events out of our sample of 95 CMEs, we could identify an associated GOES 
class flare. Figure [9] shows the distribution of the GOES classes of the associated flares, 
which were not partially occulted (61 events): 9 events are associated with an A class flare, 
27 with a B class flare, 20 events with a C class flare, and 5 events with an M class flare. 
None of the CME events under study is associated with an X class flare. 

In the upper panels of Figures [3H8] the distributions of CMEs associated with flares 
(70 events) are overplotted in color. On average, the CME peak velocity Vmax is higher 
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for events associated with flares than for those with no flare association. Assuming a 
lognormal distribution we find fi* = 495^^^^ km s~^ for events associated with flares and 
fi* = 406^^'^3 km s~^ for events with no flare-association. All CMEs, which reached a 
velocity higher than 1 000 km s~^ are associa ted with a flare. These results are in accordance 



with former studies. 



Burkepile et al. 



(120041 1 ■ who measured the average velocity of 111 



limb CMEs observed with the SMM coronagraph/polarimeter, found mean values of 566 



± 67 km s ^ 



Vrsnak et al 



'or CMEs with and 444 ± 59 km s ^ for CMEs without associated flares. 



( 120051 ) analysed the CME mean velocity and the velocity at a distance of 3 Rq 
of CMEs using LASCO data and found for the flare associated CME sample higher mean 
values than for the whole sample of CMEs. However, the Kolmogorov-Smirnov test applied 
to our sample does not show a clear distinction of both subsamples up to a significance level 
of 0.2. 

For the CME peak acceleration amax the differences in mean values and median values 
for events associated with flares and events with no flare association are even higher than 
for the velocities (cf. Table [1]). We find a more than twice as high fj,* = 551^9Q4 m s~^ for 
events with flare association than for CMEs without a flare, fi* = 2A0^1H m s~^. Although 
the range of Umax values for flare associated CMEs is large (77-6781 m s~^), all events 
which reached a peak acceleration higher than 1600 m s^^ (9 events out of our sample of 95 
CMEs) were associated with a flare. The Kolmogorov-Smirnov test shows for the parameter 
flmax at a significance level of 0.05 that both distributions (with/without associated flare) 
do not come from the same population. 

For the acceleration phase duration, we found fi* = 25.3^2ao ™^ events with 
associated flares, compared to the fi* = 4:3.8^12 4 min for events with no associated flare. All 
CMEs with tacc <14 min were associated with a flare. Applying the Kolmogorov-Smirnov 
test it can be shown that both distributions do not come from the same population at a 
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0.10 significance level. 

The measured CME height parameters Hq, h^^nax and hamax are on average smaller for 
CMEs associated with flares than for CMEs with no flare (see Table [1]). /i* of Hq and h^rnax 
are more than twice as high for events with no flare than for events with an associated 
flare. For /i* of hamax the differences between these subgroups are lower (ratio ~ 1.32). The 
Kolmogorov-Smirnov test showed for all height parameters that both subgroups do not 
come from the same population at a significance level of 0.10 for ho, 0.15 for h^max and even 
0.05 for h 

am.n.T. • 



3.1.2. CME - Filament Relation 

24 out of 95 CMEs were associated with an erupting filament. In the bottom panels of 
Figures [3H8] the distributions of these events are overplotted in color. Events, which could 
be associated with filaments showed lower v^ax than events without filament eruptions. For 
the former group we find ji* = 420^ ^1 km s~^, whereas for the latter group of CMEs we 



obtain /i* of Vmax of 531^}34 km s ^. iBurkepile et al.l (12004 ) found no significant difference 



between both subgroups for the mean velocity of ~520 km s ^. 

Similar to the CME-flare association, the differences in the mean and median values 
for Umax are more significant than for Vmax- CMEs which could be associated with filaments 
showed lower amax, with fi* of 28l;;j;264 ^ None of these events reached amax higher than 
1600 m s~^, 75% have amax <400 m s~^. CME events without associated filament eruption 
show fi* of 543;,;g"^5 m s~^, i.e. the difference between these two subgroups is about a factor 
of 2. amax was the only parameter, for which the Kolmogorov-Smirnov test suggested that 
both distributions (filament/non- filament associated events) do not come from the same 
population at a significance level of 0.05. 
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-24.00 events with 



For the CME acceleration duration tacc we found /i* of 33.1]^gyQQ 
associated and fi* of 27.8;^24 74 events with no filament association. Because filament 

eruptions are correlated to CMEs with larger a^ax it is not surprising that they are also 



correlated to events with 
and tacc (see 



Vrsnak et al 



onger 



20071 : 



tacc, due to the strong negative correlation between 



Bein et al 



201ll ). 



All CME height parameters ho, h^max and hamax have on average larger values for 
CMEs associated with filament eruptions than for CMEs without filament eruptions (Table 
d]). The differences in the median values of the lognormal distribution are about a factor of 
1.5. 



Considering the mean and median values of CME parameters we found higher values for 
Vmax and Umax and smaller values for tacc, ho, h^max and hamax for events with associated 
flares than for those with no flare. This effect is smallest for the CME peak velocity. The 
differences in the mean and median values show that CMEs associated with flares tend to 
have a more impulsive and intense acceleration and start from lower heights in the corona. 
CMEs which could be associated with an erupting filament behave oppositely. Their mean 
and median values of Vmax and a^nax are lower than for CMEs with no erupting filament 
and the mean and median values of tacc, ho, h^max and hamax are higher, when a filament 
eruption was associated. The Kolmogorov-Smirnov test shows for all CME parameters but 
for Vmax that the distributions of fiare/non-fiare associated events do not come from the 
same population. Considering filament /non-filament events the Kolmogorov-Smirnov test 
shows only for a^ax a distinctions between both distributions. 
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all 
CMEs 
(95 events) 


with 
GOES flare 
(70 events) 


without 
GOES flare 
(25 events) 


with 
filament 
(24 events) 


without 
filament 
(71 events) 


Vmax [km S'^] 


mean ± s 
median ± mad 
li±a 


526 ± 263 
460 ± 160 
6.15 ± 0.22 


550 ± 268 
461 ± 150 
6.20 ± 0.12 


459 ± 242 
393 ± 135 
6.00 ± 0.25 


467 ± 227 
400 ± 121 
6.04 ± 0.21 


546 ± 273 
461 ± 161 

6.19 ± 0.22 


a-max [m s-'^] 


mean ± s 
median ± mad 


757 ± 1034 
414 ± 246 
6.10 ± 1.05 


896 ± 1148 
556 ± 306 
6.31 ± 0.97 


367 ± 424 
241 ± 134 
5.48 ± 0.85 


391 ± 380 
266 ± 103 
5.64 ± 0.66 


880 ± 1152 
545 ± 321 
6.28 ± 1.00 


tacc [min] 


mean ± s 
median ± mad 

// ± (T 


44.6 ± 60.4 
29.0 ± 14.5 
3.28 ± 1.04 


33.9 ± 30.2 
26.0 ± 13.5 
3.23 ± 0.58 


74.4 ± 102.0 
33.0 ± 14.5 
3.77 ± 1.06 


63.1 ± 102.3 
35.0 ± 17.0 
3.50 ± 1.29 


38.3 ± 36.1 
26.0 ± 16.0 
3.33 ± 0.64 


ho [Rq] 


mean ± s 
median ± mad 

± IT 


0.24 ± 0.29 
0.14 ± 0.08 
-1.85 ± 0.88 


0.19 ± 0.24 
0.12 ± 0.06 
-2.17 ± 0.99 


0.40 ± 0.35 
0.27 ± 0.16 
-1.19 it 0.57 


0.31 ± 0.33 
0.23 ± 0.17 
-1.56 ± 0.77 


0.22 ± 0.28 
0.14 ± 0.07 
-1.97 ± 0.93 


h-amax [-^0] 


mean ± s 
median ± mad 

// ± (T 


0.53 ± 0.64 
0.26 ± 0.12 
-1.10 ± 0.90 


0.36 ± 0.37 
0.24 ± 0.10 
-1.38 ± 0.72 


0.99 ± 0.95 
0.52 ± 0.34 
-0.34 ± 0.65 


0.72 ± 0.87 
0.35 ± 0.21 
-0.79 ± 0.90 


0.46 ± 0.53 
0.26 ± 0.12 
-1.20 ± 0.84 


hvmax [-^0] 


mean ± s 
median ± mad 
± cr 


1.56 ± 1.82 
0.78 ± 0.42 
-0.09 ± 0.94 


1.17 d= 1.48 
0.70 ± 0.35 
-0.32 ± 0.96 


2.26 ± 2.38 
1.05 ± 0.53 
0.44 ± 0.75 


2.00 ± 2.42 
0.86 ± 0.40 
0.24 ± 0.91 


1.28 ± 1.54 
0.76 ± 0.41 
-0.21 ± 0.90 



Table 1: Arithmetic mean with standard deviation, median with median absolute deviation 
(mad) and mean value fi with standard deviation a of the lognormal probability function 
of the following decisive CME parameters: peak velocity Vmax, peak acceleration amax, 
acceleration phase duration tacc, height ho where the CME leading edge could be identified 
for the first time, height at peak acceleration hamax, height at peak velocity h^max- Mean, 
median and /i of all these parameters were derived for the whole sample of events (third 
column) and for several subgroups: events with associated fiares (fourth column), events for 
which no flare could be associated (flfth column), events with associated fllament eruptions 
(sixth column) , events for which no filament eruptions could be associated (seventh column) . 
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3.2. Relationship between CME and flare parameters 

We analysed the scaling of the GOES peak flux Fsxr and the GOES SXR rise time 
tsxR with the characteristics of the associated CMEs, i.e. Vmax, cimax^ tacc, ^q, h^max and 
hamax and calculated their dependency in logarithmic space. For these correlations we only 
considered those 61 (out of 70) flare events, which were not partially occulted from Earth 
view. The correlations are in general weak. Figures [10] and [11] show four representative 
correlations. 



We used a bootstrap method (lEfron]ll979l ) within a Mo nte Carlo algorithm to est imate 



19861). A 



the standard error of our calculated correlation coefficients (iEfron fc Tibshiranil 
bootstrap sample is constructed by randomly drawing with replacement from the original 
data pairs, with the bootstrap sample having the same number of entries than the original 
sample. Such bootstrap samples are constructed 1000 times and for each bootstrap sample 
the correlation coefficient is calculated. The mean value plus standard deviation of all these 
1000 calculations is used as a robust measure of the correlation coefficient and its standard 
error. 

We found a positive correlation between Fsxr and Vmax of c = 0.32±0.13 (bottom 
panel of Figure [TOjl . i.e. CMEs which were associated with flares with highe r SXR flux 



tend to reach higher peak velocities. This is consistent with the results from iMaricic et al. 



(120071 ) who analysed the correlation between SXR flares and CMEs for a set of 18 events, 
including B to X class flares, flnding c =0.52. The power law dependency between Fsxr 
and Vmax of our study can be expressed by: 

_ in3-18±0.18 p0.08±0.03 (^\ 
'Jmax — -"-U SXR ' \ J 

The correlation between Fsxr and amax is c =0.28±0.12, see the top panel of Figure [TOl i.e. 
CMEs with high peak accelerations tend to be associated with flares of high GOES SXR 
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flux. iMaricic et al.l (120071 ) also found a positve correlation between these two parameters 



with a correlation coefficient of c=0.60. We found a power law dependency of: 

n — in3-62±0.38 p0.14±0.06 ( 
'^raax — -"-U SXR " \ J 

We determined the SXR rise time tsxR from the derivative of the GOES flux curves as the 
time difference between the start of the increase and the maximum of the SXR derivative, 
which was possible for 57 events. The distribution of the time differences is shown in Figure 
[T2l The GOES SXR rise time tsxR and CME parameters (Figure [TTj) do not yield distinct 
correlations. The highest correlation coefficient was obtained for the relation between tsxR 
and the CME acceleration duration, c = 0.37±0.15. The power low dependency between 
these two parameters can be expressed as: 

Between amax and tsxR we found a negative correlation with c = —0.32 ± 0.15 and a power 
law dependency of 

„ _ in4-02±0.48 j.-0.43±0.16 (o\ 
(^max — -LU ''SXR i v"J 

i.e. CMEs with high peak accelerations are preferentially associated with flares, which show 
short SXR rise times. 

One possible explanation for the weak flare-CME correlations we obtained is the lack 
of strong flare events in our sample which contains only a few M class flares and no X class 
flares due to the solar activity minimum condition in the period under study (2007-2010). 
Since we cover only low energetic flares and thus a relatively narrow range of GOES classes, 
possible correlations may be hidden in the scatter of the data points. It turned out that 
weak flares show a larger scatter in the distribution of values which weakens our correlation 
results since predominantly A and B class flares are studied in our sample. To circumvent 
this effect we calculated the mean and median value of the CME parameters separately 
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for each GOES class (shown in Table [2]). By this procedure the CME parameters were 
weighted for each GOES class equally which reduced the influence of the scattering of the 
large number of A- and B-class flares in our sample. 

Figure [13] shows these mean values for Vmax, dmax and tacc plotted against the mean 
GOES flux of each subgroup. In this representation a clear trend is noticeable that CMEs, 



which reach higher peak velocities are related to flar es with higher GO 



Moon et al 



jS peak flux (top 



(120021 ) analysed 3217 



panel in Figure [T3|) . consistent with former studies. 
LASCO CMEs, which were associated with flares and compared the median value of the 
CME speeds derived from the whole sample, from CMEs associated with flares >C1 and 
from CMEs associated with flares >M1. When considering the whole sample the authors 
found the lowest median value, whereas CMEs, which were associated with flares >M1, 



Vrsnak et al 



(120051) 



showed the largest median value for the CME speed distribution, 
compared the velocities between B and C flare associated CMEs with the velocities of M- 
and X-class flare associated CMEs and found that the mean value of the M- and X-class 
flare associated events is about 1.4 times of the mean value of the other subgroup. 

We also flnd that the mean values of the CME peak acceleration is positively correlated 
with the GOES peak flux (middle panel of Figure [T3|) . The bottom panel of Figure [13] shows 
the mean values of the CME acceleration duration against the GOES peak flux, which 
shows a tendency that tacc is smaller for CMEs associated with higher energetic flares. The 
mean and median values calculated for each GOES class separately for all CME parameters 
under study are summarized in Table [2] 
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A-class 


B-class 


C-class 


M-class 






(9 events) 


(27 events) 


(20 events) 


(5 events) 




mean ± s 


501±258 


532±257 


609±301 


754±340 




median ± mad 


461±221 


460±135 


530±200 


710±122 




mean ± s 


597±474 


986±1236 


1229±1489 


1230±906 




median ± mad 


404±92 


639±370 


810±502 


1221±411 




mean ± s 


32.6±15.6 


28.9±24.1 


25.2±15.7 


29.9±17.7 




median ± mad 


26.0±13.0 


21.0±12.0 


20.5±9.5 


20.5±7.5 


ho 


mean ± s 


0.22±0.19 


0.19±0.27 


0.14±0.19 


0.27±0.39 




median ± mad 


0.19±0.09 


0.11±0.05 


0.10±0.05 


0.11±0.03 




mean ± s 


1.03±0.62 


0.89±1.01 


0.70±0.47 


0.80±0.65 




median ± mad 


1.30±0.61 


0.70±0.28 


0.61±0.21 


0.70±0.20 


^amax 


mean ± s 


0.41±0.28 


0.31±0.33 


0.31±0.32 


0.58±0.58 




median ± mad 


0.33±0.17 


0.23±0.08 


0.22±0.10 


0.27±0.13 



Table 2: Mean value with standard deviation and median with the median absolute deviation 
(mad) of the CME parameters Vmax, ^max, tacc, ho, Kmax and hamax, separately for the flare 
GOES classes. 
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3.3. Flare-CME temporal relationship 

We studied the relative timing of the CME acceleration phase and the flare energy 
release in order to investigate possible temporal synchronizations between these two 
phenomena. From the CME observations we derived the start (taccstart), peak [taccpeak] 
and end time (taccend) of the CME acceleration. From the GOES flare observations we 
derived the start and peak time of the 1-8 A SXR flux as well as the peak of the 1-8 A flux 
derivative. For this study we used only flare events, which were not partially occulted and 
for which a peak in the derivative of the GOES flux could be derived. This was the case for 
57 out of 70 flare-CME pairs. 

The top panels of Figure [H] show the histograms of the time differences between 
the GOES flare start times and the start of the associated CME acceleration. For the 
distribution on the left hand side the differences expressed in minutes were used, whereas 
the values on the right hand side were normalized by the CME acceleration duration. For 
the time differences we found a mean value of 4.9 ± 10.2 min and a median of 4.6 ± 6.3 
min. The distribution shows that the values are conflned to within about ±20 min (93% 
of the events lie within this interval). In 75% of the cases, the time difference is positive, 
i.e. the CME acceleration starts before the rise of the flare soft X-ray flux. Considering 



the normalized values, the time difference is sma! 



smaller than ±0.25 for 44%. 



Maricic et al. 



ler than ±0.5 for 77% of the events and 



( 120071 ) also found that for the majority of the 



events the CME acceleration starts before the SXR flare with a mean value of 23 ± 30 min 
(for a set of 18 events containing a signiflcant fraction of gradual events). This mean value 
is larger than in our study, but the time delays normalized by tacc are quite similar in bot h 



studies, with a mean value of 0.17 ± 0.71 in our study and 0.14 ± 0.20 in 



Maricic et al. 



fl2007h . 

In the middle panels of Figure [TH we compare the end of the CME impulsive 
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acceleration phase tacc_end with the GOES SXR flare peak time, which relates to the end 
of the impulsive flare energy release phase. We found that for the majority of the events 
(77%), the GOES peak occurred earlier than taccend- The mean value of the distribution 
is —9.3 ± 15.8 min and the median —6.5 ± 7.5 min. For some gradual CMEs {tacc >50 
min), the delay between flare peak and acceleration end can be >30 min (3 events). At 
the right hand side the normalized difference between the GOES peak flux and taccend is 
plotted. The distribution shows an asymmetric behaviour with a tail at positive values, and 
a distinct mode in the range of —0.5 to —0.3. For 65% of the events, the normalized time 
differences are between —0.7 and —0.1. 

The time difference between the peak of the GOES 1-8 A flux derivative and the CME 
acceleration peak (bottom left panel of Figure [T^ is found to be small. The distribution 
does strongly peak aro und zero; for the me an value we obtained 1.3 ± 8.6 min, for the 



median 1.0 ± 4.0 min. iMaricic et al.l ( l2007l ). who used the logarithmic derivative of the 
SXR GOES flux, found a mean value of 2.7 ± 14 min, consistent with our findings. In 
81% of the events, we found a difference smaller than ±10 min, for 58% of the events it is 
smaller than ±5 min. Only one event of our study showed a time delay >20 min. 65% of 
the normalized values (distribution shown at the right hand side of Figure [T^ lie within 
±0.25. The statistical parameters for the time differences described above are summarized 
in Table [31 

These findings provide strong evidence that the timing of the flare energy release and 
the CME dynam ics (acceleration) a re well synchronized, supporting previous results from 



case studies (e.g. 



Zhang et al. 



20041 ). Recent studies also analysed the relationship between 



the flare HXR emission and the CME acceleration peak, fin ding that the flare HXR pe ak 



occurs close in time with the maximum CME acceleration fiTemmer et al. 



2008 



2010 



Cheng et al. 



201G 



Berkebile-Stoiser et al. 



2OI2I). 
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Time Delay 


Minimum 


Maximum 


Arithmetic Mean ± 
Standard Deviation 


Median 
± Mad 


GOES flare start time — taccstart [min] 


-17.3 


32.6 


4.9 ± 10.2 


4.6 ± 6.3 


(GOES flare start time - taccstart) /tacc 


-1.58 


3.26 


0.17 ± 0.71 


0.17±0.20 


GOES flare peak time — tacc.end [min] 


-71.5 


21.5 


-9.3 ± 15.8 


-6.5 ± 7.5 


(GOES flare start time - tacc.end)/tacc 


-1.08 


2.83 


-0.13 ± 0.69 


-0.33 ± 0.20 


GOES der. peak time — tacc.peak [min] 


-31.5 


23.0 


1.29 ± 8.63 


1.00 ± 4.00 


(GOES der. peak time - tacc.peak) /tacc 


-0.63 


1.17 


0.09 ± 0.38 


0.02 ± 0.18 



Table 3: Statistical parameters (minimum value, maximum value, arithmetic mean, standard 
deviation, median and median absolute deviation) for time delays between characteristic 
CME and flare times: time difference between the start of the GOES flare and the start 
of the CME acceleration taccstart (differences expresssed in minutes in the first row, values 
normalized by the CME acceleration time tacc in the second row), time delay between the 
GOES flare peak time and the end of the CME acceleration phase taccend (delays expressed 
in minutes in the third row, normalized values in the forth row) and time difference between 
the peak time of the GOES derivative and the peak time of the CME acceleration taccpeak 
(delays expressed in minutes in the fifth, normalized values in the sixth row). 
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4. Summciry and Conclusion 

Based on a sample of 95 CME events we present a statistical study on various 
characteristic CME parameters and their relation to flares. CMEs which are associated 
with flares show on average higher peak velocities {vmax), higher peak accelerations {amax), 
shorter acceleration phase durations (tacc), lower heights at peak velocity (hymax), lower 
heights at peak acceleration (hamax) and lower initiation heights (ho). The ratio between 
the median values of the lognormal probability functions of both subgroups is about a 
factor of 2. Only for Vmax the ratio is signiflcantly smaller (~1.1), most probably due to 
the small range of Vmax values. Due to the anticorrelation between amax and tacc and the 
relation Vmax ~ o,max • tacc the range for Vmax is smaller than for the other two quantities, 
it basically covers only one order of magnitude. Although we found clear differences in the 
mean and median values of the two subgroups there exist also events associated with flares, 
which have low v^ax and Umax values and high tacc, h^max, hamax and ho. For instance, the 
smallest measured amax value for a CME with flare was 77 m s^^ and the highest measured 
value for a CME without flare was 1577 m s~^. The Kolmogorov-Smirnov test suggested a 
distinction between both distributions (flare/non-flare associated events) for every CME 
parameter (except for Vmax) at 0.05-0.15 levels of significance. The clearest distinctions 
were found for amax and hamax at a 0.05 level of significance. 

CMEs which are associated with erupting filaments show on average smaller Vmax, 
ttmax, and larger tacc, h^max, hamax and ho- These trends are the other way round than for 
the flare association. The ratio between the mean and median values of both subgroups is 
somewhat smaller than for the flare association. Again we found the smallest ratio between 
the median values //* of Vmax (l-l)- The fi* of amax and h^max showed the highest ratio with 
1.9. For tacc, hamax and ho we found ratios between 1.2 and 1.5. But there exist also CMEs 
associated with fllament eruptions, which have high Vmax and amax values and low tacc. 
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hvmax, hamax and hf). FoT example the highest measured amax value for a CME which was 
associated with an erupting filament was 1561 m s~^, the smallest value for a CME event 
with no erupting filament was 35 m s~^ (the second lowest value of the whole distribution). 
Both of these events were not associated with a flare. The Kolmogorov-Smirnov test did 
not show a clear distinction between both subgroups. Only for amax the test suggested that 
both distributions do not come from the same population at a 0.05 level of significance. 

The correlations obtained between Vmax, (^max and tacc with the GOES peak flux 
F'sxR and the SXR rise time tsxR of the associated flare were low. We found a weak 
positive correlation between the CME acceleration duration tacc and the flare rise time tsxR 
(c = 0.37±0.15) and a weak negative correlation between amax and tsxR (c = — 0.32±0.15). 
Correlations between Vmax and amax with Fsxr. were c = 0.32 ± 0.13 and c = 0.28 ± 0.12 
respectively. If the events are averaged and binned into the different GOES classes, 
the correlations are much more distinct. The mean values in each GOES class show an 
increasing trend for Vmax and amax, and a decreasing trend for tacc with higher GOES flux. 

For the majority of the events (75%) we found that th e CME acceleration starts before 



Maricic et al. 



(120071 ). suggesting 



the SXR flare onset, which is consistent with the findings of [ 
that the flare is a consequence of the eruption. Similar to our study these authors also 
found about one fourth of the events, for which the onset of the associated flare occurred 
before the CME acceleration started. They explained these cases by a superposition of 
two flares, a conflned flare in the pre-eruption stage, which releases only a part of the 
stored magnetic energy, and a second flare, beginning after the CME acceleration onset 
and causing a prolongation of the flrst flare in the full-disk integrated SXR light curve. The 
CME is associated with the second flare but because of the superposition, the flare start is 
measured from the flrst one. To test this hypothesis, we checked the SXR curves for all 
events, for which the CME acceleration start was after the flare onset and found indeed in 
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I 



11 out of 14 events evidence of a second SXR peak, confirming their suggestion Q|. Figure 
[TS] shows two representative examples. The top panels show the CME acceleration profile, 
the bottom panels the GOES SXR flux together with its derivative. We marked the two 
possible SXR peaks by arrows. Assuming that we have a superposition of two subsequent 
flares, we would also measure erroneous SXR start times and as a result too long tsxR- This 
misidentification would also influence our correlation negatively. To test this, we correlated 
tsxR again with a^ax and tacc considering only events for which the CME acceleration starts 
before the flare onset and found indeed distinctly higher correlation coefficients (Figure [16], 
c = 0.59 ± 0.12 and c = —0.50 ± 0.14) than for the whole sample (shown in Figure [TT| 
c=0.37±0.15 and c = —0.32 ± 0.15). Thus superposition of two subsequent flares and more 
complex structures are probably a reason for the weak correlations and may account for a 
considerable number of events, where the flare seems to start before the eruption. 

For the majority of the events (77%) we found that the end of the CME acceleration 
occurred after the SXR peak. Especially long duration flares reach a certain point, when 
they become too weak to compensate cooling of the hot plasma. As a result the SXR curve 
decreases although the energy release in the flare may still be going on. 

For 81% of the events the time delay between the CME acceleration peak and the 
peak of the GOES SXR flux derivative, which is a proxy for the flare energy release rate. 



2004 



Zhang fc Dere 



Temmer et al. 



2006 : 



Maricic et al. 



2007 



CME and the flare energy release ( 


Lin 




Temmer et al. 


2008; 


Reeves & Moats 


2010 



2010l ). There are basically two forces acting on a flux rope in equilibrium. 



an upward directed magnetic pressure and a downward directed magnetic tension of 



^This is a significantly higher rate than in the total sample of events, in which about 20% 
showed a double SXR peak 
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the overlying magnetic field. When the magnetic structure looses equilibrium, it starts 



rising and a current s. 



( jPriest fc Forbes 



leet is formed below it, where magnetic reconnection takes place 



20021 ). The reconnection reduces the tension of the overlying field and 



enhances the magnetic pressure at the bottom part of the fl ux rope due to additional 



20081 ). The upward 



poloidal flux, providing the upward acceleration of the rope ( IVrsnakl 
motion of the rope leads to elongation of the current sheet and a more efficient reconnection, 
thus enhancing the acceleration. On the other hand, more efficient reconnection means 
also a more powerful energy release in the CME-associated flare, which directly relates the 
dynamics of the eruption and the energy release in the flare. 
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Fig. 1.— CME kinematics and GOES 1-8 A soft X-ray flux for the CME-flare events tfiat 
occurred on 2008 January 7 (left) and 2009 January 9 (right). The top panels show the 
measured CME height-time curve derived from STEREO EUVI (crosses), CORl (triangles) 
and C0R2 (asterisks) observations. The measurement errors (0.03 Rq for EUVI, 0.125 Rq 
for CORl and 0.3 Rq for C0R2), which are in some cases smaller than the plot symbols, and 
the spline fit (solid line) are overplotted. The second and third panels show the CME velocity 
and acceleration profiles, derived from numerical differentiation of the direct measurements 
(plot symbols) and the spline fit (solid line) to the height-time curve. The grey shaded area 
represents the error range of the spline fit. The bottom panels show the GOES flux (black 
solid line) and its derivative (red dashed line) of the associated flare. 
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Fig. 2. — Same as Figure [T]but for the events observed on 2009 December 22 (left) and 2010 
February 12 (right). 
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Fig. 3. — Histogram of the CME peak velocity v^ax for the whole sample of 95 events (grey 
distributions in the top and bottom panels). In the top panel the Vmax distribution of CME 
events associated with a flare is overlaid in color, in the bottom panel the histogram of CME 
events associated with an erupting filament is overlaid. 
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Fig. 4. — Same as Fig. |3]but for the CME peak acceleration. 
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Fig. 5. — Same as Fig. |3]but for the CME acceleration duration. 
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Fig. 6. — Same as Fig. [3] but for the height Hq where the CME could be first identified. This 
height is used as an estimate of the CME initiation height. 
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Fig. 7. — Same as Fig. |3]but for the height h^max-, where the CMEs reached their maximum 
velocity. 
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Fig. 8. — Same as Fig. [3] but for the height hamax, where the CMEs reached their maximum 
acceleration. 



30 \ - 

25 : : 

20 - 1 - 

15 : : 

10 - : 

5 : 1 : 

\ I - 

A B C M X 

Fig. 9. — Distribution of the GOES classes of the flare events under study. 
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Fig. 10. — CME peak velocity (top) and peak acceleration (bottom) against the GOES 
peak flux of the associated flare. A double logarithmic space is used for the plot and the 
calculation of the correlation coefficient c, which is annotated in each panel. The regression 
line is overplotted in red. 
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Fig. 11. — CME acceleration duration (top) and peak acceleration (bottom) against the flare 
SXR rise time. The correlation coeflicient c (calculated in double logarithmic space) and the 
regression line is overplotted. 
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Fig. 12. — Distribution of the risetime tsxR of the GOES SXR flares. 



-34- 



^ 800 


X 




'</> 


c = 0.90 




E 


X 




X 


X 




o 
E 


X 




^ 400 






10"^ 10"' 10"^ 10"^ 10"^ 


2000 








c = 0.95 X X 




^ 1000 








X 




X 

o 


X 




E 

o 






300 






10"^ 10"' 10"^ 10"^ 10""* 


35 








X 




c 


X 










O 


X 




U 

o 


X 




-1—' 


c = -0.86 




28 






10"^ 10"' 10"^ 10"^ 10""* 




GOES peok flux [W m"'] 




Mean value of the CME peak velocity v^ax (top), acceleration 



acceleration duration tacc (bottom) for each GOES class plotted against the mean GOES flux 
in each subgroup. The correlation coefficient c (calculated in logarithmic space) is annotated 
in each panel. 
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Fig. 14. — Distribution of the time differences between the start of the GOES flare and 
the start of the CME acceleration (top), GOES peak time and CME acceleration end time 
(middle) and the peak time of the derivative of the GOES 1-8 A flux and the CME accelera- 
tion peak time (bottom). On the left hand side, the time differences are plotted in minutes, 
whereas on the right hand side the time differences are normalized by the acceleration du- 
ration of the corresponding CMEs. Positive values mean that the CME acceleration start, 
end or peak occurred before the flare SXR start, peak or derivative peak, respectively. 
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Fig. 15. — CME acceleration (top) and GOES SXR curve (bottom) from the events observed 
on 25 October 2009 (left hand side) and 23 March 2008 (right hand side). In both events 
the CME acceleration starts after the flare onset, which may be related to an overlap of two 
flares. Arrows mark two SXR maxima indicating two possible flare maxima. 
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Fig. 16. — Same as Figure [TT] but only for events, where the CME acceleration starts before 
the flare. 
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